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ABSTRACT 

We study the expansion of the nearby Universe using a sample of Type la 
supernovae at redshifts below 0.08. Both the conventional supernova standard- 
ization method weighing heavily on the magnitude at optical maximum, and 
the color based CMAGIC method are employed to put the supernovae on the 
Hubble diagram. These supernovae allow peculiar velocities of nearby galaxies 
to be measured at unprecedented precision. We have investigated in detail the 
possibility of a varying Hubble constant with redshift and found no evidence 
of a monopole term for the nearby Universe. A large scale streaming motion 
is found at an amplitude of about 340 63 71 km/sec, aligned in the direction of 
(7o, bo) = (312°. 0l_ 3 7 5 4 , 25°. 7^g 2 ), which is close to the direction of the center of 
Shapley supercluster of galaxies. This streaming motion is best fit by a function 
involving a strong bipolar term. The streaming velocity field extends from the 
lowest redshift (~ 0.007) to beyond 0.025 and likely out to even higher redshifts. 
The velocity field at redshift below 0.01 can be equally well described by a dipole 
field or by the same bipolar streaming velocity field that reaches out to beyond 
z ~ 0.025. We are also able to deduce a robust estimate of the random velocity 
component of the peculiar velocity field. Within the volume of redshift below 
0.01 (weighted average redshift of ~ 0.067), this thermal component is found to 
be about 270 km/sec. After correcting this smooth streaming motion, we are 
able to significantly improve the Hubble expansion fits of these supernovae. The 
CMAGIC method gives a dramatic decrease of x 2 from 90 to 63 for 69 degrees of 
freedom, and yields a residual scatter of only 0.12 magnitude; the correction also 
gives a moederate improvement to Hubble fit for the maximum light method. A 
whole sky Type la supernova search with carefully controlled systematic error is 
the key to confirm the large scale velocity patterns identified in this paper. The 
large scale streaming motion of the local Universe is of fundamental importance 
as it probes quantitatively the statistical isotropy of the local Universe. The 
bipolar pattern identified in this paper, if born out with future observations, will 
shed light on the formation of the structures of the Universe from the smallest 
to the largest scale. 

Subject headings: Supernovae - Cosmology - Distance scales 
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Introduction 



Type la supernovae (SNIa) form a uniform class of o bjects th a t are preci se cosmic dis- 
tance indicators. The conventional method, established by iPhillipsI (119931 ) and lHamuy et al 
( 119961 ). is based on the empirical correlation between the brightness at optica l maximum 
of SN Ia and the decline rate of the light curve after optical maximum (see also iRiess et al. 
( 119961 ); iPerlmutter et al.l (119991 )). Multi-color light curves of SNIa are required to correct 
for the extinction caused by dust on the line of sight to the supernovae in the host galax- 
ies. Dust extinction correction uncertainty usually dominants the error budget, not only 
because of the technical difficulties in acquiring multi-color light curves, but also because 
of the lack of knowledges of the extinction properties of dust particles in extragalactic en- 
vironment. Even with the simplified assumption that the extinction law derived for dust 
particles in the Galaxy also applies to extragalactic environment, the coefficients of absolute 
extinction to reddening are known to vary from galaxies to galaxies. For example, the value 
R B = A(B)/E(B — V) has been deduced to be in a range of as low as slightly below 2 
for some galaxies - a value which is significantly lower than the ca nonical value of 4.3 as 
has been found for th e average of interstellar dust in the Galaxy (jElias-Rosa et al.l 12006 ; 
Krisciunas et al.ll2007l ). Furthermore, extinction correction may be complicated by the pres- 
ence of circu mstellar dust particles in which scattered light may affect the observed light and 
color curves (jWangll2005l ). Because these methods rely heavily on the magnitude at optical 
maximum, we will refer to this conventional method as maximum light method (MLM). 

The Color-Magnitude Intercept Calibration (CMAGIC) of SNIa p rovides an alternativ e 
approach that is l ess sensitive to dust extinction I Wang et al.l (l2003bl ); IConley et al.l (120061 ); 
Wang et al.l (120061 ) . In CMAGIC, the magnitude at a given color index is used as the primary 
distance indicator. The empirical CMAGIC linear relation can be employed to deduce the 
magnitude at a given color index. The CMAGIC linear relation shows that on the color- 
magnitude diagram, from about one week to one month past optical maximum (the exact 
epoch s of the linear relation is a function of the light curve width, see IWang et al.l (J2003b, 
20061 )). the color- magnitude diagram of SNIa can be well approximated by a straight line 
with a universal slope. For example, for B and V band data, this relation is given as 
B = Bq — (3(B — V), with f3 = 1.94. This simple relation reduces the degree of freedom 
of the light curve fits, and when compared to the MLM method, yields significantly reduced 
statistical errors (jWang et al.l 120061 ; IConley et al.l 120061 ) in the final distance modulus. In 



Conley et al.l (120061 ). application of CMAGIC to high redshift SNIa shows that with a small 
sub-set of SNe with enough data for the application of CMAGIC, it is possible to achieve 
constraints on cosmological parameters that are tighter than those derived from a bigger 
sample of SNe but using the conventional method. 



The velocity field of the local universe has been explored previously by several groups 
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(IZehavi et al.l Il998l ; lHaugboelle et al.l 120061 ) using distance mo duli of SNIa. Th e se studies 



are b ased mostly on the MCLS approach of the MLM method (iRiess et a 
2006 ). A tentativ e mon opole field is identified by IZehavi et al.l (119981 ) and 
Haueboelle et al.l (120061 ) recently modeled the distance moduli published in 



1998 



Jha et al 



ha et al 



(2007) 



Jha et all (12006 



2007h and identified dipole and quadruple components of the local peculiar velocity field. 



Large scale streaming patterns of the local universe can be a source of systematic errors 
for supernova cosmology if not well understood. A streaming motion of about 300 km/sec can 
introduce a magnitude bias larger than 0.02 mag out to redshift of 0.1. Since all high redshift 
SNIa surveys need to anchor the Hubble diagram to a low redshift sample, a streaming 
peculiar velocity field can introduce a noticeable bias for future precision cosmology studies 
targeting accuracies of the order of 0.01 mag. 

In this study, we apply the MLM and CMAGIC methods to a sample of nearby SNIa 
to study the expansion of the local universe at z < 0.08. Our goal is to study the angular 
distribution of the peculiar velocities of the host galaxies to probe the kinematics of the local 
universe. In this paper, §2] describes the data sample and the CMAGIC Hubble diagram, §3] 
analyzes the streaming motion and decomposes them into monopole, dipole, and quadruple 
components, §H gives the Hubble diagram after corrections of the streaming velocity flow, 
and a brief discussion of the current findings is given in §51 



2. The Data Sample and the CMAGIC Hubble diagrams 



The detail s of the SN sample , the light curve fit and the CMAGIC method are described 
i n Table 2 of (Wang et al. 2006h. but this s tudy adds a f ew additional SNe: SN 2003du 



( Anupama et al 



(IKrisciunas et al 



200 



5 



Stanishev et al. 



20071), SN 2005am (ILi et al 



2007 



SN 2004eo (IPastorello et al.ll2007ah. 2004 S 



20061 ). and SN 2005cf (IPastorello et al.ll2007M ) 



As in IWang et al.l (120061 ). we have used the the SuperStretch algorithm to fit the super- 
nova light curves. A major advantage of both SuperStretch and CMAGIC methods is that 
the light or color curve models can be applied to the raw observed data before extinction 
and filter K corrections. Unlike other light curve fitting methods which assume that all 
SNIa light curves can be described by a one parameter family of functions, the SuperStretch 
method is a multi-paramter light curve fitting algorithm with a maximum of 6 free parame- 
ters of each light curve of which 4 are relevant to the light curve sh apes. The SuperStretch 
method is built upon t he success of the conventional stretch method (IPerlmutter et al.lll999 ; 
Goldhaber et al.ll200ll ). but allows the stretch parameter to vary with time from optical max- 
imum, and introduces a secondary component to compensate the light curve bumps in R and 
J, and sometime bumps in B and V bands for some peculiar SNe. As such, most of the light 
curve parameters such as magnitudes at B maximum, Amis, and CMAGIC parameters are 
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identical to those published in IWang et al.l (120061 ). As shown in IWang et al.l (l2003bl . 120061 ). 
very frequently the CMAGIC linear fits give x 2 /DoF that are smaller than 1, suggesting 
that the statistical noise of the original data is often very conservative, or with unknown 
correlations. In this study, the errors of light curve and CMAGIC fits are scaled to \ 2 P er 
degree of freedom to be one, whereas a slightly more sophisticated approach is adopted in 



Wang et al.l ( 2006T). This only introducs some small differences between the errors in this 



paper and in lWang et al.l (120061 ). Such differences are not expected to significantly alter the 
statistical properties of the data set. 

This paper studies only SN Ia with suffic i ent ob servations to allow for our multi-parameter 
light curve and CMAGIC fits (jWang et al. ( 2006 )). In addition, the following criteria are 
applied to select the final data sample: (a) Only SNIa with Am i 5 smaller than 1.7 are cho- 
sen - this practically excludes all SN 1991bg-like fast decliners ( IFilippenko et al.lll992l ); (b) 
SNIa with B max — V max > 0.5 are excluded - this reduces potential errors in extinction 
correction of heavily reddened SNe; and (c) only SNe at redshifts below 0.1 are selected. 
The original data set has a total of 125 SNIa of which 76 satisfy the above criterions. The 
redshift distribution of the sample is shown in Figure HJ Most of the SNe are at redshift 
below 0.05. 

The magnitudes deduced from SuperStretch and CMAGIC fits need to be corrected by 
the light curve shapes. We will analyze in parallel both the MLM and CMAGIC method in 
this paper. For the MLM in B and V bands, the distance modulus to a supernova is given 
by 



B n 



Mi - ai(Ara 15 - 1.1) - bi (B r , 



\ 7( 



where M.\ is the absolute B magnitude at optical maximum, a\ and b 2 are correction coeffi- 
cients to be determined together with parameters of cosmological interest, B max , V max , and 
Am 15 are the magnitude at B and V maximum, and the 5-band magnitude decline during 
the first 15 days past B maximum, respectively. The CMAGIC method does not make use 
of the magnitude at V max . As shown in an our earlier studies, it is equivalent to 



At = B BV - M 2 - a 2 {Am 15 - 1.1) - (b 2 - P)[ 



B r , 



^ + , + 1 ,(I-^),, (2) 



where M. 2 is the absolute B BV magnitude, a 2 and b 2 are coefficients that needs to be de- 
termined simultaneously with other cosmological parameters, Bbv and /3 are the CMAGIC 
intercept magnitude and slope (IWang et al.ll2003bl ). the number 0.6 inside the square bracket 
is because B BV is calculated at B — V = 0.6. As shown in previous studies, bi and b 2 
has a range of value s from 1.5-4 depending on details of the subsets of supernovae but is 
around 2 in average JTripplll997l : kstier et alJbood : IWang et al.ll2006l ). and j3 = 1.94. It is 
readily seen from equation (T5]) that the CMAGIC method is essentially independent of the 
measured colors of the SNe. 
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Fig. 1. — Histogram of the redshift distribution of the SNe included in this study. 



In this study we also deduce the R B values for Galactic extinction of SNe based on a 
library of well observed SN spectra. The R B value of SNIa evolves with time t after explosion 
due to the time evolution of SN spectra. This leads to a small correction to the Amis values of 
reddened SNe, of the order of (R B (0)-R B (15)) xE(B-V), where R B (0) and R B (15) are the 
absolute to rel ative extinction ratio at day and day 15 for Galactic dust. For the standard 
extinction law (ICardeUi et al.lll989l ) with R B = 4.1, typical values are: R B (0) = 4.15±0.06, 
R v (0) = 3.14±0.02,i2 B (15) = 4.05±0.11,andi2 v (0) = 3.12±0.16. Note that in general 
Rb — Rv 7^ 1- We also deduce estimates of the errors of R B and Ry based on the various 
spectra of SNe we have in the library. Typical values of R B and Ry uncertainties are 0.06 
and 0.02 at optical maximum, and are 0.11 and 0.16 at day 15, respectively, Although not a 
dominant source of errors, the uncertainties in the extinction coefficients are propagated to 
the errors of the distances we deduced from these SNe. The second term on the right hand 
side of equations (1) and (2) corrects the extinction due to extragalactic dust or any intrinsic 
ma gnitude dependence on the measured color. This approach, first proposed in iTrippI (119971 ) 
and lTripp fc Branch! (119991 ) . has also been employed in recent high redshift supernova studies 
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Astier et al.l (120061 ). The method has the obvious shortcomings of not being able to correctly 
separate potential systematic effects of the magnitude- intrinsic color dependence from the 
errors of dust extinction. However, this is still the best one can do with SNIa considering 
the limitations on the current understandings of the physical properties of SNIa and the 
properties of extragalactic dust. Other approaches for handling dust extinction include: (a) 
assumption of a correlation between intrinsic supernova color and the light c urve shape, 
and a theoretical prio r on the distribution of the extragalactic dust reddening (IRiess et al. 



relation as shown in 



19981 ; I Jha et al.l 12007); (b) assump t ion t hat the late time color of SNIa follows a linear 



maximum (IWang et al. 



Phill ips et al.l (119991 ); (c) using the SN color at 12 days past optical 



2005T); and (d) using near-IR observations to pinpoin t the reddening 



and dust extinction law (IKrisciunas et al.l 12004 . 120071 : lElias-Rosa et al.l 120061 ) . Not all there 
methods are applicable to the current data set. Some of these methods introduce additional 
assumptions that are difficult to justify. 

A random thermal peculiar velocity component of 350 km/sec, assumed to be Gaussian, 
is added to each supernova. This is again a poorly known quantity. However, the assumed 
peculiar velocity is comp atible with pairwise velocity dis t ribution functions de termined from 
galaxy redshift surveys (IDavis et al.l Il997l ; iLandyl 120021 ; lHawkins et al.l 120031 ) , and is likely 
to be a good guess of the thermal component of the galaxy peculiar velocity. This will be 
discussed further in the next sections. 

The quantities that are relevant to this study are shown in Tables [T] and [21 where the 
columns are: (1) the identification of the supernova, (2) the redshift in the rest frame of 
the CMB, (3) and (4) the B— band maximum magnitude (Tabled]), or the color magnitude 
(Table [2]), and error, (5) and (6) the Ami^^B) value derived using SuperStretch fit and error, 
respectively, (7) and (8) the color term as in Equations ([1]) or (T5]), and error, respectively, 
(9) the deduced distance modulus, (10) and (11) the residual on the Hubble diagram, and 
error, respectively, (12) and (13) the deduced peculiar velocity and error, respectively, (14) 
and (15) show the R.A. and Deck of the supernovae in ecliptic coordinates. For column (7) 

This 



and S(B -V) = 



5^ + 0.6 + 1.2(1-^. 



and (8), the color index is B max — V m 

definition of EiB — V) is different from that in IWang et al.l (l2003bl ) by an extra last term 
which helps to reduce the covariance between EiB — V) and (3. It is important to note that 
the quantities such as magnitude at maximum, Am 15 , and CMAGIC magnitude and color 
shown in Tables [1] and [2] are in general correlated. It is impossible to tabulate the covariance 
matrix here but they are included in the analyses throughout this paper. 

By minimizing the magnitude scatter on the Hubble diagram, we can set constraints on 
the coefficients a\, b\, a 2 , and b 2 of Equations ([!]) and (jSJ. The results are shown in Figure [2j 
We found that b 2 is arou nd 2 which is very c lose to being equal to the mean slope < (3 > of 
the CMAGIC linear fits ( IWang et al.ll2003bl ). and a 2 around 0.26. With these values, B B v 
is an excellent standard candle requiring only a very small Amis and color correction. The 
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errors listed in columns (11) an d (13) of Table [2] are in general a f e w tim es smaller than 



those given in the gold sample of iRiess et al 



and in 



Jha et al.l (120071 ). This is mainly 



due to the introduction of the CMAGIC method which reduces significantly the light curve 
shape and color dependence of the derived distance scales. 

It is frequently assumed that SNIa are not perfect standard candl es but with some in- 
trinsic magnitude dispersion ( Perlmutter et al. 1999 ; Riess et al~1ll998l ). The exact level of 
intrinsic dispersion of SNIa magnitude is of fundamental interest in applying SNIa to cosmo- 
logical studies. Earlier studies have shown that the total (intrinsic-l- m e asurement) dispersio n 
deduced for the sample of nearby supernovae (e.g. IRiess et all Jl996h : IPhillips et aP Jl999h ) 
is about 0.18 mag. The fact that the y 2 per degree of fre edom is generally we ll around 1 on 



the Hubble diagram in the fits of IRiess et al.l (119961 ) and IPhillips et al.l (Il999l ) suggests that 
0.18 mag is only an upper limit of the intrinsic dispersion. The actual intrinsic dispersion 
should be much smaller than 0.18 mag as it is highly possible that observational or dust 
extinction errors are the dominant sources of errors in those ana l yses. The intrinsic magni- 
tude dispersion of SNIa has been studied in detail in IWang et al.l (120061 ) where the observed 
magnitude scatter around the best fit Hubble diagram is at tributed to the com bination of 
a random component and a coherent correlated component. IWang et al.l (120061 ) found that 
after systematic corrections are performed, MLM and CMAGIC exhibit mutual rms intrin- 
sic variation equal to 0.074 ± 0.019 mag, of which at least an equal share likely belongs to 
CMAGIC. In this study, we probe the streaming velocity field of the local universe which 
is likely to be a dominant source of correlated noise on the Hubble diagram of the nearby 
Universe. 




Fig. 2. — (a, left): The correction coefficients a± and b\ of Equation ([TJ, and (b, right): The 
correction coefficients a 2 and b 2 for Equation ([2]). The contours are 1, 2, and 3-cx levels from 
the inner most to the outer most. This shows that errors of the light curve shape correction 
coefficients are anti-correlated for both the MLM and CMAGIC method. 



Table 1. MLM Parameters of the Supernovae 



SN 


z cm6 


Mag 


a Amis 


a Color a 


Mo 


An 


a Vpec (km/sec) 


a (km/dec) RA (2000.0) 


Dec (2000.0) 


bl 


(1) 


(2) 


(3) 


(4) (5) 


(6) (7) (8) 


(9) 


(10) 


(11) (12) 


(13) (14) 


(15) 


(16) 



1989B 


0.00359 


12.23 


0.02 


1 


.25 


0.02 


0.35 


0.03 


30.39 


-0.52 


0.06 


229.72 


23.19 


11 


20 


13.91 


+13 


00 


19.20 


31.59 


1990af 


0.04992 


17.77 


0.01 


1 


.59 


0.02 


0.04 


0.02 


36.56 


-0.14 


0.04 


944.88 


250.71 


21 


34 


58.12 


-62 


44 


7.38 


19.88 


1990N 


0.00447 


12.54 


0.03 


1 


.13 


0.02 


-0.09 


0.03 


31.94 


0.55 


0.06 


-386.85 


48.37 


12 


42 


56.70 


+13 


15 


23.69 


38.95 


1990O 


0.03060 


16.22 


0.04 





.9.3 


0.03 


-0.03 


0.04 


35.57 


-0.04 


0.09 


160.87 


378.63 


17 


15 


35.92 


+16 


19 


25.80 


15.43 


1990T 


0.03967 


17.26 


0.09 


1 


.24 


0.03 


-0.14 


0.13 


36.72 


0.53 


0.31 


-3317.41 


2166.39 


19 


59 


2.28 


-56 


15 


30.00 


26.21 


1990Y 


0.03870 


17.63 


0.31 


1 


.15 


0.05 


0.17 


0.32 


36.32 


0.19 


0.57 


-1061.54 


3382.79 


03 


37 


22.64 


-33 


02 


40.09 


-14.60 


1991ag 


0.01388 


14.45 


0.05 





.89 


0.03 


0.03 


0.09 


33.67 


-0.20 


0.21 


363.48 


368.60 


20 


00 


8.66 


-55 


22 


3.38 


25.67 


1991S 


0.05606 


17.76 


0.24 


1 


.30 


0.05 


-0.22 


0.26 


37.42 


0.45 


0.50 


-3880.42 


4797.62 


10 


29 


27.79 


+22 


00 


46.41 


17.36 


1991U 


0.03243 


16.49 


0.10 


1 


.17 


0.04 


-0.09 


0.11 


35.87 


0.13 


0.23 


-617.58 


1082.72 


13 


23 


22.20 


-26 


06 


28.70 


79.49 


1992A 


0.00594 


12.62 


0.02 


1 


.31 


0.01 


0.11 


0.02 


31.39 


-0.61 


0.06 


438.21 


34.56 


03 


36 


27.40 


-34 


57 


31.59 


-12.95 


1992ae 


0.07461 


18.53 


0.08 


1 


.24 


0.07 


0.12 


0.09 


37.30 


-0.31 


0.19 


3011.66 


1653.42 


21 


28 


17.66 


-61 


33 


0.00 


19.66 


1992ag 


0.02700 


16.25 


0.05 


1 


.07 


0.05 


0.10 


0.06 


35.19 


-0.15 


0.14 


523.87 


478.32 


13 


24 


10.12 


-23 


52 


39.30 


77.31 


1992al 


0.01350 


14.48 


0.02 


1 


.21 


0.02 


-0.05 


0.02 


33.73 


-0.07 


0.05 


129.53 


94.19 


20 


45 


56.45 


-51 


23 


40.00 


17.79 


1992au 


0.06035 


18.35 


0.55 


1 


.44 


0.07 


0.19 


0.61 


36.82 


-0.31 


1.14 


2434.49 


8609.03 


00 


10 


40.27 


-49 


56 


43.31 


-1.75 


1992bc 


0.01960 


15.11 


0.02 





.75 


0.01 


-0.03 


0.02 


34.58 


-0.04 


0.05 


109.09 


130.38 


03 


05 


17.28 


-39 


33 


39.69 


-11.02 


1992bg 


0.03648 


16.75 


0.04 


1 


.14 


0.03 


-0.02 


0.05 


35.96 


-0.04 


0.12 


216.53 


571.38 


07 


41 


56.53 


-62 


31 


8.81 


32.89 


1992bh 


0.04509 


17.61 


0.03 





.97 


0.06 


0.07 


0.03 


36.69 


0.21 


0.07 


-1384.05 


489.30 


04 


59 


27.55 


-58 


49 


44.19 


14.55 


1992bk 


0.05885 


18.26 


0.06 


1 


.63 


0.03 


0.05 


0.07 


37.02 


-0.06 


0.15 


475.81 


1187.83 


03 


43 


1.90 


-53 


37 


56.81 


4.43 


1992bl 


0.04223 


17.41 


0.06 


1 


.50 


0.04 


0.03 


0.07 


36.27 


-0.06 


0.16 


345.25 


897.23 


23 


15 


13.25 


-44 


44 


34.50 


-3.10 


1992bo 


0.01723 


15.75 


0.02 


1 


.63 


0.02 


-0.01 


0.03 


34.64 


0.30 


0.06 


-778.01 


166.24 


01 


21 


58.44 


-34 


12 


43.50 


-19.29 


1992J 


0.04612 


17.83 


0.11 


1 


.53 


0.12 


0.10 


0.16 


36.50 


-0.03 


0.35 


159.88 


2189.69 


10 


09 


0.30 


-26 


38 


24.41 


46.14 


1992P 


0.02649 


16.07 


0.02 


1 


.25 


0.04 


-0.02 


0.02 


35.22 


-0.07 


0.06 


252.86 


216.55 


12 


42 


48.95 


+10 


21 


37.50 


41.77 


1993ag 


0.05004 


17.80 


0.03 


1 


.27 


0.03 


0.08 


0.04 


36.67 


-0.04 


0.08 


263.66 


544.50 


10 


03 


35.00 


-35 


27 


47.59 


48.32 


1993ah 


0.02850 


16.41 


0.14 


1 


.46 


0.04 


-0.21 


0.19 


35.94 


0.49 


0.43 


-2149.05 


2144.05 


23 


51 


50.27 


-27 


57 


47.00 


-21.40 


1993B 


0.07008 


18.44 


0.06 


1 


.31 


0.04 


0.01 


0.07 


37.48 


0.01 


0.15 


-79.17 


1428.12 


10 


34 


51.38 


-34 


26 


30.00 


54.16 


1993H 


0.02513 


16.71 


0.02 


1 


.56 


0.02 


0.18 


0.03 


35.14 


-0.03 


0.06 


107.52 


218.93 


13 


52 


50.34 


-30 


42 


23.31 


82.81 


19930 


0.05293 


17.62 


0.03 


1 


.27 


0.03 


-0.06 


0.03 


36.85 


0.01 


0.08 


-98.95 


564.51 


13 


31 


7.88 


-33 


12 


50.50 


86.75 


1994D 


0.00261 


11.76 


0.01 


1 


.57 


0.01 


-0.10 


0.02 


30.93 


0.71 


0.04 


-301.80 


22.37 


12 


34 


2.37 


+07 


42 


4.70 


43.79 


1994M 


0.02431 


16.32 


0.03 


1 


.41 


0.03 


0.06 


0.04 


35.17 


0.07 


0.09 


-244.09 


307.86 


12 


31 


8.61 


+00 


36 


19.90 


50.33 


1994Q 


0.02987 


16.52 


0.08 





.96 


0.04 


0.09 


0.09 


35.54 


-0.02 


0.20 


69.60 


830.10 


16 


49 


51.06 


+40 


25 


55.53 


0.73 


1994S 


0.01611 


14.79 


0.02 





.88 


0.07 


-0.01 


0.04 


34.12 


-0.07 


0.10 


146.34 


219.04 


12 


31 


21.82 


+29 


08 


4.19 


22.88 


1995ak 


0.02198 


16.06 


0.05 


1 


.48 


0.05 


0.00 


0.07 


35.02 


0.14 


0.16 


-443.22 


533.96 


02 


45 


48.83 


+03 


13 


50.10 


-52.69 


1995al 


0.00588 


13.32 


0.02 





.96 


0.01 


0.13 


0.02 


32.22 


0.24 


0.05 


-206.55 


49.02 


09 


50 


55.97 


+33 


33 


9.41 


3.05 


1995bd 


0.01443 


15.22 


0.02 





.73 


0.03 


0.31 


0.05 


33.79 


-0.16 


0.09 


312.09 


166.42 


04 


45 


21.24 


+11 


04 


2.50 


-39.72 


1995D 


0.00766 


13.23 


0.02 





.97 


0.01 


-0.03 


0.02 


32.58 


0.02 


0.05 


-16.29 


50.52 


09 


40 


54.79 


+05 


08 


26.60 


21.94 



Table 1 — Continued 



SN 


z cm£> 


Mag 


a Amis 


a Color a 


Mo 




a Vp ec (km/sec) 


a (km/dec) RA (2000.0) 


Dec (2000.0) 


bl 


(1) 


(2) 


(3) 


(4) (5) 


(6) (7) (8) 


(9) 


(10) 


(11) (12) 


(13) (14) 


(15) 


(16) 



1996M 


0.03485 


16.68 


0.03 





.84 


0.04 


0.05 


0.04 


35.87 


-0.02 


0.09 


113.36 


419.58 


00 


36 


17.97 


+11 


23 


40.50 


-61.96 


1996bo 


0.01632 


15.84 


0.03 


1 


16 


0.03 


0.32 


0.03 


34.12 


-0.10 


0.07 


221.28 


145.48 


01 


48 


22.86 


+11 


31 


15.10 


-64.80 


1996C 


0.03007 


16.62 


0.08 





.94 


0.04 


0.09 


0.09 


35.65 


0.08 


0.19 


-354.78 


816.04 


13 


50 


48.59 


+49 


19 


7.09 


4.13 


1996X 


0.00800 


12.99 


0.01 


1 


.28 


0.01 


-0.02 


0.01 


32.13 


-0.53 


0.02 


519.91 


19.76 


13 


18 


1.12 


-26 


50 


45.31 


79.94 


1996Z 


0.00974 


14.41 


0.05 


1 


.05 


0.04 


0.42 


0.08 


32.50 


-0.59 


0.18 


694.31 


185.91 


09 


36 


44.82 


-21 


08 


51.70 


37.05 


1997bp 


0.00944 


13.91 


0.01 


1 


.19 


0.01 


0.15 


0.01 


32.64 


-0.38 


0.02 


458.71 


21.38 


12 


46 


53.75 


-11 


38 


33.20 


63.18 


1997bq 


0.00960 


14.44 


0.03 


1 


.16 


0.02 


0.10 


0.04 


33.32 


0.27 


0.09 


-374.69 


131.13 


10 


17 


5.33 


+73 


23 


2.12 


-24.69 


1997cw 


0.01595 


15.86 


0.07 





.99 


0.03 


0.28 


0.11 


34.37 


0.20 


0.25 


-451.71 


611.57 


00 


25 


17.27 


+12 


53 


6.20 


-62.06 


1997E 


0.01333 


15.11 


0.02 


1 


.38 


0.03 


0.05 


0.03 


33.99 


0.21 


0.06 


-408.07 


123.10 


06 


47 


38.10 


+74 


29 


51.00 


-37.88 


1998aq 


0.00426 


12.31 


0.01 


1 


.17 


0.01 


-0.11 


0.01 


31.75 


0.47 


0.02 


-305.78 


13.42 


11 


56 


25.87 


+55 


07 


43.19 


-3.85 


1998bu 


0.00416 


12.13 


0.01 


1 


.01 


0.01 


0.35 


0.01 


30.42 


-0.81 


0.02 


389.92 


9.15 


10 


46 


46.01 


+ 11 


50 


7.50 


27.98 


1998dk 


0.01204 


14.84 


0.32 


1 


.30 


0.07 


0.07 


0.35 


33.72 


0.17 


0.67 


-287.18 


1224.73 


00 


14 


32.16 


00 44 10.90 


-48.68 


1998ec 


0.02012 


16.31 


0.05 


1 


.05 


0.03 


0.17 


0.07 


35.05 


0.37 


0.17 


-1113.86 


553.99 


06 


53 


6.11 


+50 


02 


22.09 


-32.88 


1998es 


0.00957 


13.84 


0.01 





.74 


0.01 


0.08 


0.01 


33.02 


-0.03 


0.03 


37.49 


36.04 


01 


37 


17.60 


+05 


52 


50.60 


-59.40 


1998V 


0.01717 


15.10 


0.03 


1. 


.10 


0.02 


-0.03 


0.04 


34.35 


0.02 


0.07 


-44.96 


178.95 


18 


22 


37.41 


+15 


42 


8.41 


4.04 


1999aa 


0.01525 


14.74 


0.02 





.81 


0.02 


0.02 


0.02 


34.04 


-0.03 


0.04 


63.09 


91.66 


08 


27 


42.03 


+21 


29 


14.80 


-2.17 


1999ac 


0.00982 


14.09 


0.01 


1 


.23 


0.03 


0.01 


0.01 


33.16 


0.06 


0.03 


-80.55 


41.33 


16 


07 


15.02 


+07 


58 


20.40 


32.58 


1999aw 


0.03924 


16.69 


0.02 





.78 


0.02 


-0.05 


0.02 


36.19 


0.03 


0.04 


-159.23 


244.86 


11 


01 


36.37 


-06 


06 


31.60 


44.16 


1999dk 


0.01395 


14.82 


0.03 


1 


.00 


0.03 


0.03 


0.03 


33.97 


0.09 


0.07 


-183.78 


148.41 


01 


31 


26.92 


+14 


17 


5.70 


-67.83 


1999dq 


0.01308 


14.41 


0.02 





.94 


0.03 


0.06 


0.02 


33.51 


-0.23 


0.05 


388.27 


85.65 


02 


33 


59.68 


+20 


58 


30.41 


-69.48 


1999ee 


0.01055 


14.85 


0.01 





.90 


0.01 


0.29 


0.01 


33.38 


0.12 


0.03 


-172.99 


40.38 


22 


16 


10.00 


-36 


50 


39.69 


-3.79 


1999ek 


0.01760 


15.57 


0.01 


1 


.08 


0.01 


0.11 


0.06 


34.48 


0.09 


0.08 


-223.52 


214.10 


05 


36 


31.60 


+16 


38 


17.80 


-32.86 


1999gh 


0.00879 


14.51 


0.07 


1 


.67 


0.03 


0.31 


0.07 


32.53 


-0.33 


0.14 


371.78 


149.36 


09 


44 


19.75 


-21 


16 


25.00 


38.64 


2000ca 


0.02452 


15.58 


0.01 





.93 


0.01 


-0.03 


0.01 


34.93 


-0.19 


0.03 


613.89 


85.22 


13 


35 


22.98 


-34 


09 


37.00 


87.61 


2000cf 


0.03603 


17.09 


0.05 


1 


.31 


0.04 


-0.02 


0.05 


36.19 


0.21 


0.10 


-1124.49 


560.47 


15 


52 


56.19 


+65 


56 


12.69 


-15.94 


2000cn 


0.02321 


16.60 


0.01 


1 


.58 


0.02 


0.20 


0.01 


34.97 


-0.03 


0.03 


86.30 


103.97 


17 


57 


40.42 


+27 


49 


58.09 


0.43 


2000dk 


0.01645 


15.35 


0.02 


1. 


.43 


0.03 


0.02 


0.03 


34.28 


0.04 


0.06 


-100.00 


128.02 


01 


07 


23.52 


+32 


24 


23.19 


-83.49 


2000E 


0.00422 


12.78 


0.01 


1 


.03 


0.02 


0.03 


0.04 


31.93 


0.67 


0.07 


-455.01 


58.20 


20 


37 


13.77 


+66 


05 


50.19 


-39.39 


2000fa 


0.02179 


15.75 


0.05 


1 


.09 


0.02 


-0.10 


0.06 


35.20 


0.34 


0.12 


-1125.79 


412.65 


07 


15 


29.88 


+23 


25 


42.41 


-16.84 


2001ba 


0.03053 


16.19 


0.01 


1 


.06 


0.01 


-0.08 


0.01 


35.60 


-0.00 


0.03 


15.85 


144.11 


11 


38 


1.76 


-32 


19 


51.00 


66.20 


2001bt 


0.01445 


15.27 


0.01 


1 


.29 


0.01 


0.17 


0.01 


33.89 


-0.06 


0.02 


116.43 


38.54 


19 


13 


46.75 


-59 


17 


22.81 


32.94 


2001el 


0.00365 


12.79 


0.00 


1 


.16 


0.00 


0.12 


0.01 


31.62 


0.68 


0.01 


-400.54 


7.76 


03 


44 


30.57 


-44 


38 


23.69 


-3.55 


2001V 


0.01604 


14.61 


0.03 





.73 


0.03 


0.06 


0.07 


33.84 


-0.34 


0.18 


701.75 


332.88 


11 


57 


24.93 


+25 


12 


9.00 


24.42 


2002bo 


0.00529 


13.94 


0.01 


1. 


.24 


0.01 


0.42 


0.01 


31.91 


0.15 


0.02 


-117.46 


13.09 


10 


18 


6.51 


+21 


49 


41.70 


15.89 


2002cl 


0.02238 


16.15 


0.03 


1 


.37 


0.02 


-0.03 


0.03 


35.25 


0.33 


0.07 


-1116.58 


265.68 


20 


56 


30.09 


-18 


33 


34.30 


-4.89 



Table 1 — Continued 



SN 


z cmf> 


Mag 




Ami 5 


a 


Color 


a 


Mo 


A/i 




Vp ec (km/sec) 


a (km/dec) 


RA (2000.0) 


Dec (2000.0) 


bf 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


(14) 


(15) 


(16) 


2002er 


0.00855 


14.26 


0.01 


1.26 


0.01 


0.15 


0.02 


32.96 


0.15 


0.04 


-186.42 


49.59 


17 11 29.88 


+07 59 44.80 


22.08 


2003du 


0.00665 


13.50 


0.04 


1.15 


0.03 


-0.09 


0.05 


32.91 


0.65 


0.11 


-698.57 


132.64 


14 34 35.80 


+59 20 3.81 


-6.66 


2004eo 


0.01473 


15.08 


0.00 


1.33 


0.00 


0.07 


0.01 


33.94 


-0.05 


0.02 


106.82 


36.45 


20 32 54.19 


+09 55 42.70 


-18.07 


2004S 


0.00986 


14.24 


0.02 


1.18 


0.01 


0.08 


0.02 


33.16 


0.05 


0.05 


-66.36 


71.79 


06 45 43.50 


-31 13 52.50 


9.78 


2005am 


0.00897 


13.69 


0.01 


1.54 


0.02 


0.09 


0.02 


32.38 


-0.53 


0.04 


584.00 


43.53 


09 16 12.47 


-16 18 16.00 


30.36 


2005cf 


0.00704 


13.20 


0.00 


1.14 


0.00 


-0.01 


0.01 


32.39 


0.01 


0.02 


-11.42 


16.69 


15 21 32.20 


-07 24 47.50 


51.73 



Table 2. CMAGIC Parameters of the Supernovae 



SN 


z cm6 


Mag 


a Amis 


a Color a 


Mo 


An 


a Vp ec (km/sec) 


a (km/dec) RA (2000.0) 


Dec (2000.0) 


bl 


(1) 


(2) 


(3) 


(4) (5) 


(6) (7) (8) 


(9) 


(10) 


(11) (12) 


(13) (14) 


(15) 


(16) 



1989B 


0.00359 


12.59 


0.02 


1 


.25 


0.02 


0.39 


0.02 


30.48 


-0.43 


0.03 


192.42 


10.81 


11 


20 


13.91 


+13 


00 


19.20 


31.59 


1990af 


0.04992 


18.98 


0.03 


1 


.59 


0.02 


-0.01 


0.03 


36.69 


-0.01 


0.03 


75.54 


230.29 


21 


34 


58.12 


-62 


44 


7.38 


19.88 


1990N 


0.00447 


13.83 


0.01 


1 


.13 


0.02 


-0.05 


0.02 


31.67 


0.28 


0.02 


-185.38 


14.51 


12 


42 


56.70 


+13 


15 


23.69 


38.95 


1990O 


0.03060 


17.57 


0.01 





.93 


0.03 


-0.07 


0.02 


35.46 


-0.14 


0.06 


589.60 


227.05 


17 


15 


35.92 


+16 


19 


25.80 


15.43 


1990T 


0.03967 


18.38 


0.00 


1 


.24 


0.03 


0.04 


0.04 


36.19 


0.01 


0.03 


-34.41 


185.27 


19 


59 


2.28 


-56 


15 


30.00 


26.21 


1990Y 


0.03870 


18.22 


0.06 


1 


.15 


0.05 


0.29 


0.15 


36.11 


-0.02 


0.04 


132.32 


235.90 


03 


37 


22.64 


-33 


02 


40.09 


-14.60 


1991ag 


0.01388 


15.80 


0.05 





.89 


0.03 


-0.08 


0.05 


33.72 


-0.14 


0.07 


268.57 


122.44 


20 


00 


8.66 


-55 


22 


3.38 


25.67 


1991S 


0.05606 


19.07 


0.02 


1 


.30 


0.05 


-0.08 


0.17 


36.90 


-0.06 


0.08 


489.33 


628.09 


10 


29 


27.79 


+22 


00 


46.41 


17.36 


1991U 


0.03243 


17.64 


0.02 


1 


.17 


0.04 


0.02 


0.05 


35.48 


-0.26 


0.04 


1092.84 


170.85 


13 


23 


22.20 


-26 


06 


28.70 


79.49 


1992A 


0.00594 


13.65 


0.02 


1 


.31 


0.01 


0.08 


0.03 


31.45 


-0.56 


0.03 


402.61 


15.97 


03 


36 


27.40 


-34 


57 


31.59 


-12.95 


1992ae 


0.07461 


19.81 


0.05 


1 


.24 


0.07 


-0.06 


0.10 


37.65 


0.03 


0.07 


-335.62 


739.98 


21 


28 


17.66 


-61 


33 


0.00 


19.66 


1992ag 


0.02700 


17.22 


0.10 


1 


.07 


0.05 


0.12 


0.08 


35.07 


-0.26 


0.11 


903.74 


377.14 


13 


24 


10.12 


-23 


52 


39.30 


77.31 


1992al 


0.01350 


15.92 


0.01 


1 


.21 


0.02 


-0.11 


0.02 


33.72 


-0.09 


0.02 


157.58 


40.70 


20 


45 


56.45 


-51 


23 


40.00 


17.79 


1992au 


0.06035 


19.45 


0.05 


1 


.44 


0.07 


0.06 


0.36 


37.19 


0.06 


0.12 


-483.97 


1036.97 


00 


10 


40.27 


-49 


56 


43.31 


-1.75 


1992bc 


0.01960 


16.81 


0.04 





.75 


0.01 


-0.25 


0.05 


34.74 


0.12 


0.06 


-336.03 


159.89 


03 


05 


17.28 


-39 


33 


39.69 


-11.02 


1992bg 


0.03648 


18.09 


0.02 


1 


.14 


0.03 


-0.06 


0.03 


35.91 


-0.08 


0.11 


419.01 


547.94 


07 


41 


56.53 


-62 


31 


8.81 


32.89 


1992bh 


0.04509 


18.73 


0.02 





.97 


0.06 


0.04 


0.02 


36.63 


0.16 


0.03 


-1001.76 


176.62 


04 


59 


27.55 


-58 


49 


44.19 


14.55 


1992bk 


0.05885 


19.38 


0.01 


1 


.63 


0.03 


0.05 


0.04 


37.06 


-0.01 


0.02 


84.65 


176.44 


03 


43 


1.90 


-53 


37 


56.81 


4.43 


1992bl 


0.04223 


18.66 


0.05 


1 


.50 


0.04 


-0.02 


0.05 


36.39 


0.06 


0.04 


-367.83 


235.79 


23 


15 


13.25 


-44 


44 


34.50 


-3.10 


1992bo 


0.01723 


16.92 


0.02 


1 


.63 


0.02 


0.01 


0.03 


34.62 


0.28 


0.02 


-719.47 


66.34 


01 


21 


58.44 


-34 


12 


43.50 


-19.29 


1992J 


0.04612 


18.67 


0.11 


1 


.53 


0.12 


0.18 


0.10 


36.41 


-0.11 


0.09 


710.27 


524.77 


10 


09 


0.30 


-26 


38 


24.41 


46.14 


1992P 


0.02649 


17.83 


0.43 


1 


.25 


0.04 


-0.21 


0.27 


35.48 


0.19 


0.15 


-746.38 


596.00 


12 


42 


48.95 


+10 


21 


37.50 


41.77 


1993ag 


0.05004 


18.86 


0.03 


1 


.27 


0.03 


0.07 


0.04 


36.66 


-0.04 


0.07 


303.88 


504.07 


10 


03 


35.00 


-35 


27 


47.59 


48.32 


1993ah 


0.02850 


17.70 


0.03 


1 


.46 


0.04 


-0.06 


0.10 


35.47 


0.02 


0.05 


-75.12 


211.84 


23 


51 


50.27 


-27 


57 


47.00 


-21.40 


1993B 


0.07008 


19.61 


0.05 


1 


.31 


0.04 


0.01 


0.07 


37.40 


-0.07 


0.07 


695.94 


612.54 


10 


34 


51.38 


-34 


26 


30.00 


54.16 


1993H 


0.02513 


17.40 


0.02 


1 


.56 


0.02 


0.25 


0.02 


35.14 


-0.03 


0.04 


101.81 


141.08 


13 


52 


50.34 


-30 


42 


23.31 


82.81 


19930 


0.05293 


18.96 


0.03 


1 


.27 


0.03 


-0.07 


0.06 


36.75 


-0.08 


0.05 


582.89 


325.53 


13 


31 


7.88 


-33 


12 


50.50 


86.75 


1994D 


0.00261 


13.22 


0.05 


1 


.57 


0.01 


-0.13 


0.05 


30.92 


0.70 


0.04 


-299.14 


18.15 


12 


34 


2.37 


+07 


42 


4.70 


43.79 


1994M 


0.02431 


17.50 


0.17 


1 


.41 


0.03 


0.01 


0.22 


35.26 


0.17 


0.14 


-577.43 


513.17 


12 


31 


8.61 


+00 


36 


19.90 


50.33 


1994Q 


0.02987 


17.60 


0.03 





.96 


0.04 


0.06 


0.05 


35.50 


-0.05 


0.03 


223.04 


114.42 


16 


49 


51.06 


+40 


25 


55.53 


0.73 


1994S 


0.01611 


16.21 


0.01 





.88 


0.07 


-0.11 


0.02 


34.12 


-0.07 


0.03 


147.84 


61.47 


12 


31 


21.82 


+29 


08 


4.19 


22.88 


1995ak 


0.02198 


17.18 


0.02 


1 


.48 


0.05 


0.04 


0.04 


34.92 


0.05 


0.04 


-148.92 


111.91 


02 


45 


48.83 


+03 


13 


50.10 


-52.69 


1995al 


0.00588 


14.32 


0.04 





.96 


0.01 


0.10 


0.04 


32.22 


0.24 


0.04 


-203.81 


36.20 


09 


50 


55.97 


+33 


33 


9.41 


3.05 


1995bd 


0.01443 


16.15 


0.05 





.73 


0.03 


0.12 


0.06 


34.14 


0.19 


0.30 


-396.86 


656.67 


04 


45 


21.24 


+11 


04 


2.50 


-39.72 


1995D 


0.00766 


14.61 


0.02 





.97 


0.01 


-0.08 


0.03 


32.48 


-0.08 


0.04 


83.27 


39.81 


09 


40 


54.79 


+05 


08 


26.60 


21.94 



Table 2 — Continued 



SN 


z cm£> 


Mag 


a Amis 


a Color a 


Mo 




a Vp ec (km/sec) 


a (km/dec) RA (2000.0) 


Dec (2000.0) 


bl 


(1) 


(2) 


(3) 


(4) (5) 


(6) (7) (8) 


(9) 


(10) 


(11) (12) 


(13) (14) 


(15) 


(16) 



1996M 


0.03485 


17.93 


0.02 





.84 


0.04 


-0.03 


0.03 


35.87 


-0.03 


0.07 


147.37 


320.88 


00 


36 


17.97 


+11 


23 


40.50 


-61.96 


1996bo 


0.01632 


16.34 


0.02 


1 


16 


0.03 


0.33 


0.02 


34.23 


0.01 


0.05 


-19.10 


118.05 


01 


48 


22.86 


+11 


31 


15.10 


-64.80 


1996C 


0.03007 


17.85 


0.04 





.94 


0.04 


-0.02 


0.05 


35.75 


0.18 


0.04 


-797.53 


178.15 


13 


50 


48.59 


+49 


19 


7.09 


4.13 


1996X 


0.00800 


14.26 


0.02 


1 


.28 


0.01 


-0.04 


0.03 


32.06 


-0.60 


0.05 


579.17 


39.65 


13 


18 


1.12 


-26 


50 


45.31 


79.94 


1996Z 


0.00974 


14.90 


0.01 


1 


.05 


0.04 


0.28 


0.03 


32.92 


-0.17 


0.06 


218.22 


72.94 


09 


36 


44.82 


-21 


08 


51.70 


37.05 


1997bp 


0.00944 


14.64 


0.04 


1 


.19 


0.01 


0.25 


0.02 


32.45 


-0.57 


0.04 


656.95 


38.63 


12 


46 


53.75 


-11 


38 


33.20 


63.18 


1997bq 


0.00960 


15.28 


0.02 


1 


.16 


0.02 


0.18 


0.02 


33.12 


0.07 


0.02 


-87.73 


29.02 


10 


17 


5.33 


+73 


23 


2.12 


-24.69 


1997cw 


0.01595 


16.49 


0.03 





.99 


0.03 


0.31 


0.05 


34.33 


0.16 


0.05 


-377.28 


120.63 


00 


25 


17.27 


+12 


53 


6.20 


-62.06 


1997E 


0.01333 


16.23 


0.01 


1. 


.38 


0.03 


0.04 


0.03 


34.00 


0.22 


0.08 


-428.68 


157.00 


06 


47 


38.10 


+74 


29 


51.00 


-37.88 


1998aq 


0.00426 


13.85 


0.03 


1 


.17 


0.01 


-0.16 


0.03 


31.66 


0.37 


0.03 


-239.64 


20.33 


11 


56 


25.87 


+55 


07 


43.19 


-3.85 


1998bu 


0.00416 


12.85 


0.02 


1 


.01 


0.01 


0.24 


0.02 


30.74 


-0.50 


0.02 


254.75 


10.41 


10 


46 


46.01 


+ 11 


50 


7.50 


27.98 


1998dk 


0.01204 


15.86 


0.00 


1 


.30 


0.07 


0.11 


0.20 


33.61 


0.06 


0.07 


-106.50 


112.57 


00 


14 


32.16 


00 44 10.90 


-48.68 


1998ec 


0.02012 


17.08 


0.02 


1 


.05 


0.03 


0.21 


0.03 


34.97 


0.29 


0.05 


-874.46 


173.10 


06 


53 


6.11 


+50 


02 


22.09 


-32.88 


1998es 


0.00957 


15.10 


0.03 





.74 


0.01 


-0.04 


0.09 


33.07 


0.02 


0.04 


-30.06 


52.77 


01 


37 


17.60 


+05 


52 


50.60 


-59.40 


1998V 


0.01717 


16.65 


0.25 


1 


.10 


0.02 


-0.12 


0.30 


34.39 


0.06 


0.14 


-139.19 


335.96 


18 


22 


37.41 


+15 


42 


8.41 


4.04 


1999aa 


0.01525 


16.16 


0.01 





.81 


0.02 


-0.13 


0.02 


34.12 


0.05 


0.03 


-115.88 


61.86 


08 


27 


42.03 


+21 


29 


14.80 


-2.17 


1999ac 


0.00982 


15.05 


0.00 


1 


.23 


0.03 


0.09 


0.01 


32.92 


-0.19 


0.03 


241.03 


37.28 


16 


07 


15.02 


+07 


58 


20.40 


32.58 


1999aw 


0.03924 


18.24 


0.02 





.78 


0.02 


-0.20 


0.05 


36.24 


0.07 


0.05 


-404.16 


269.47 


11 


01 


36.37 


-06 


06 


31.60 


44.16 


1999dk 


0.01395 


15.88 


0.03 


1 


.00 


0.03 


0.06 


0.03 


33.78 


-0.10 


0.04 


184.53 


79.48 


01 


31 


26.92 


+14 


17 


5.70 


-67.83 


1999dq 


0.01308 


15.47 


0.09 





.94 


0.03 


0.06 


0.06 


33.38 


-0.35 


0.10 


589.50 


157.41 


02 


33 


59.68 


+20 


58 


30.41 


-69.48 


1999ee 


0.01055 


15.65 


0.02 





.90 


0.01 


0.19 


0.01 


33.59 


0.33 


0.02 


-513.26 


32.23 


22 


16 


10.00 


-36 


50 


39.69 


-3.79 


1999ek 


0.01760 


16.58 


0.04 


1. 


.08 


0.01 


0.09 


0.06 


34.45 


0.07 


0.34 


-161.92 


854.34 


05 


36 


31.60 


+16 


38 


17.80 


-32.86 


1999gh 


0.00879 


14.92 


0.02 


1 


.67 


0.03 


0.39 


0.04 


32.64 


-0.23 


0.04 


263.87 


41.83 


09 


44 


19.75 


-21 


16 


25.00 


38.64 


2000ca 


0.02452 


17.28 


0.05 





.93 


0.01 


-0.20 


0.05 


35.07 


-0.05 


0.04 


167.06 


145.52 


13 


35 


22.98 


-34 


09 


37.00 


87.61 


2000cf 


0.03603 


18.33 


0.03 


1 


.31 


0.04 


-0.02 


0.04 


36.12 


0.14 


0.04 


-737.69 


192.47 


15 


52 


56.19 


+65 


56 


12.69 


-15.94 


2000cn 


0.02321 


17.32 


0.04 


1 


.58 


0.02 


0.24 


0.09 


35.04 


0.04 


0.12 


-143.94 


378.10 


17 


57 


40.42 


+27 


49 


58.09 


0.43 


2000dk 


0.01645 


16.42 


0.04 


1 


.43 


0.03 


0.07 


0.06 


34.18 


-0.06 


0.06 


129.33 


131.96 


01 


07 


23.52 


+32 


24 


23.19 


-83.49 


2000E 


0.00422 


13.88 


0.00 


1 


.03 


0.02 


0.06 


0.04 


31.74 


0.47 


0.22 


-309.42 


160.15 


20 


37 


13.77 


+66 


05 


50.19 


-39.39 


2000fa 


0.02179 


17.28 


0.02 


1 


.09 


0.02 


-0.12 


0.03 


35.04 


0.18 


0.05 


-565.02 


148.52 


07 


15 


29.88 


+23 


25 


42.41 


-16.84 


2001ba 


0.03053 


17.66 


0.03 


1 


.06 


0.01 


-0.13 


0.05 


35.51 


-0.10 


0.05 


398.44 


218.97 


11 


38 


1.76 


-32 


19 


51.00 


66.20 


2001bt 


0.01445 


16.01 


0.02 


1 


.29 


0.01 


0.23 


0.02 


33.83 


-0.13 


0.04 


243.21 


80.00 


19 


13 


46.75 


-59 


17 


22.81 


32.94 


2001el 


0.00365 


13.62 


0.01 


1. 


.16 


0.00 


0.18 


0.01 


31.48 


0.54 


0.01 


-308.31 


6.27 


03 


44 


30.57 


-44 


38 


23.69 


-3.55 


2001V 


0.01604 


15.77 


0.06 





.73 


0.03 


0.02 


0.08 


33.73 


-0.45 


0.07 


906.79 


126.76 


11 


57 


24.93 


+25 


12 


9.00 


24.42 


2002bo 


0.00529 


14.41 


0.08 


1 


.24 


0.01 


0.40 


0.05 


32.18 


0.42 


0.05 


-340.00 


47.96 


10 


18 


6.51 


+21 


49 


41.70 


15.89 


2002cl 


0.02238 


17.38 


0.02 


1 


.37 


0.02 


-0.01 


0.03 


35.15 


0.24 


0.05 


-778.70 


187.50 


20 


56 


30.09 


-18 


33 


34.30 


-4.89 
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SN 


z cmf> 


Mag 




Ami 5 


a 


Color 


a 


Mo 


An 




Vp ec (km/sec) 


a (km/dec) 


RA (2000.0) 


Dec (2000.0) 


bf 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


(14) 


(15) 


(16) 


2002er 


0.00855 


15.06 


0.04 


1.26 


0.01 


0.20 


0.03 


32.87 


0.07 


0.10 


-84.75 


120.45 


17 11 29.88 


+07 59 44.80 


22.08 


2003du 


0.00665 


14.97 


0.02 


1.15 


0.03 


-0.13 


0.03 


32.80 


0.55 


0.02 


-575.72 


26.64 


14 34 35.80 


+59 20 3.81 


-6.66 


2004eo 


0.01473 


16.06 


0.02 


1.33 


0.00 


0.10 


0.02 


33.86 


-0.14 


0.07 


268.80 


128.20 


20 32 54.19 


+09 55 42.70 


-18.07 


2004S 


0.00986 


15.27 


0.01 


1.18 


0.01 


0.08 


0.02 


33.12 


0.00 


0.06 


-1.81 


84.18 


06 45 43.50 


-31 13 52.50 


9.78 


2005am 


0.00897 


14.62 


0.01 


1.54 


0.02 


0.14 


0.01 


32.35 


-0.55 


0.03 


606.13 


32.94 


09 16 12.47 


-16 18 16.00 


30.36 


2005cf 


0.00704 


14.46 


0.01 


1.14 


0.00 


-0.03 


0.02 


32.30 


-0.08 


0.06 


75.45 


56.78 


15 21 32.20 


-07 24 47.50 


51.73 
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Fig. 3. — Hubble diagrams from the conventional MLM method (a, upper panel) and from 
the CMAGIC (b, lower panel). Notice that the error bars for the CMAGIC are significantly 
smaller than for MLM. The CMAGIC method also yields a significantly smaller magnitude 
dispersion. The solid lines show the best fit Hubble diagram. The dotted lines show the 
magnitude deviation introduced by a random peculiar velocity of 350 km/sec. The inserts 
show the histogram of the magnitude deviations from the best fit Hubble diagram. The solid 
lines in the inserts show the ideogram of the residual distribution. The baseline cosmology 
model assumes Qm = 0.3 and fl\ = 0.7. 
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Figure [3] (a) and (b) shows the resulting Hubble diagrams for the MLM and CMAGIC 
methods, based on Equations ([I]) and (j2J), respectively. The solid lines are from cosmological 
models of a flat universe with Qm — 0.3 and Qm — 0.7. A Hubble constant of 71 
km/sec/Mpc is assumed in these fits, but the exact value of that does not affect the goodness 
of the fits. It is seen that the CMAGIC Hubble diagram has a magnitude dispersion of 0.14 
mag which is smaller than the 0.18 mag for the MLM method. This indicates that the 
CMAGIC method does help tightening up the dispersions on the Hubble diagram. However, 
the CMAGIC method yields a significantly worse Hubble fit, with a x 2 = 90.68 for 73 
degrees of freedom, which is remarkably larger than the \ 2 = 78.65 for the same degrees of 
freedom but with MLM. The probability of get x 2 larger than 90.68 from a random sample 
with 73 degrees of freedom is about 8%, suggesting that either the CMAGIC method has 
underestimated the errors or unidentified systematic effects such as extinction errors, or large 
scale streaming motions may be important. If so, such effects have shown up strongly with 
CMAGIC but the MLM method may be less sensitive to those effects. 

With the assumed peculiar velocity of 350 km/sec, the Hubble diagram of the MLM 
method allows very little room for a significant intrinsic dispersion. Adding a small intrinsic 
dispersion to the CMAGIC method will reduce the x 2 to an acceptable level, however, we 
will see below that this is unlikely to be the case; the failure of a satisfactory Hubble fit with 
CMAGIC actually carries important information about the streaming motion of the local 
universe. 



3. The Peculiar Velocity Field 

There is little physical motivations of the functional form of the streaming velocity 
field. Conventionally it has been modeled in terms of spherical harmonic expansions and the 
monopole, dipole, and quadruple terms are studied. 



3.1. The Monopole Flow 



It has been suggested that the nearby Universe may have an expansion velocity that is 
different from the Hubble expansion at higher redshift. Su ch a configuration may arise if the 
Galaxy is are inside a low density void. iKim et al.l (119971 ) compared supernova distances at 
redshifts arou nd 0.4 and 0.1 and se t a 2-cr upper limit of 10% for an outflow within the local 
300 h" 1 Mpc. IZehavi et all (119981 ). based on a larger sample of well measured supernovae, 
studied the redshift dependence of the Hubble constant, and found a putative void extending 
to about 7 h~* Mpc. The existence of such a void is confirmed recently at a 2-cr level by 
Jha et al.l (120071 ) based on an even larger sample of SNe. 



- 16 - 



To investigate the nature of the monopole field, we performed a joint Hubble diagram 
fit using the following model for the expansion of the Universe 



V L (H,q,z ,z,Q M ,Q A ) = [1/q + (1 - l/q)H(z - z )]D L (H, z, M , A ), (3) 

where Dl is the luminosity distance in the Friedmann- Robertson- Walker Universe, H the 
Hubble constant at redshift above zq, 7i(z) the Heaviside function, zq the redshift cut above 
which the flow is consistent with the Hubble flow with Hubble constant H, Qm, and ar e 
the density parameters for matter and dark energy, respectively. We assume Qm = 0.3, 
and Qa = 0.7 in these studies. 

The observed supernova magnitude is then fit by minimizing the quantity 



X 



^K(^) - 5logV L (H, q, z , z u 0.3, 0.7) - 25] 2 /af 



(4) 



where Zi is the observed redshift of the supernovae, m\ is the observed SN magnitude 
corrected by Equations ([I]) or ([2]), the errors of the corrected magnitudes. The three 
additional free parameters are H, q, and zq and the free parameters are determined simul- 
taneously through x 2 minimization. In such a fit, the parameter H actually represents the 
combined effect of the Hubble constant and the zero point of the absolute magnitude of 
the supernovae. The final constraint on Zq and q is derived by marginalizing over nuisance 
parameters such as H, and the light curve shape and color correction coefficients. 



First, we tried to repro duce the monopo le field shown in |jha et al.l (120071 ). We used 
the distance moduli given in iJha et al.l (120071 ) and applied the \ 2 minimization procedure 
described above. The sample of th e SNe was chosen to match exactly the 95 Hubble flow 
[zc > 2,500km/ sec) SNe used by IJha et al.l (120071 ). The resulting constraint on z and q 
is shown in Figure HI Indeed, at a level slightly above 2-cr the fit sug gest q = 1.046, at 
zq = 0.026, which is in agreement with the Hubble bubble reported in I Jha et al.l (120071 ). 



The presence of a Hubble bubble at such a large distance scale is of sufficient cosmological 
significance that it deserves careful examination using different approaches. As is well known, 
the uncertainties in extinction correction is by far the dominant source of errors in supernova 
cosmology. To further inve stigate the authe nticity of the Hubble bubble, we have analyzed 
the extinction correction of I Jha et al.l (120071 ). Figure [5] shows t he magnitude res iduals on the 



Hubble diagram and the extinction in ^/-band (Ay) as given in I Jha et al.l (120071 ). A linear fit 
to these data points reveals a significant correlation, as shown in Figure [5j The \ 2 reduced 
from 117.3 to 101 after the introduction of the linear fit which suggests a 87% confidence 
level. The deduced relation is given by Residual = 0.077(±0.025) + 0.336(±0.085)A y , 
where the numbers in parenthesis are the errors of the corresponding coefficients. This 
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Fig. 4. — (a) L eft hand p a nel, c onfidence level on the monopole field using the distance 
moduli given in Ijha et al.l (120071 ). (b) Right hand panel, the confidence level calculated 
using the same data as in (a) but after correcting the linear relation shown in Figure [51 The 
horizontal axes are the cut-off redshift, and the vertical axes the ratio of Hubble constant 
below and above z minus 1, as defined in Equation (J3j). The contours show the 1-a (solid 
lines), 2-<T (dotted lines), and 3-a (dashed lines) significance levels. 



residual correlation with Ay suggests that the extinction correction in Ijha et al.l (120071 ) is 
systematically incorrect for SNe with different amount of reddening. Such a correlation 
violates the requirement that the distance to an SN be independent of the reddening to the 
SN. 

As there is not yet a reliable model to determine the color-excess of SNIa, the alternative 
approach is to completely remove the correlation of the distance and any non-distance related 
parameters such as Amis and the measured color. This is the approach used in equations 
(CE]) and (j2|) . As an exercise based on the same principle, we may correct the distance moduli 
in (jjha et al.l 120071 ) with the linear correlation shown in Figure [5J The resulting distance 
moduli are again fitted by the broken Hubble law given in Equation (J3j). The confidence 
levels of the ratio of the Hubble constants below and above the cut-off redshift zq are shown 
in Figure H](b). This shows a dramatic change as compared to Figure H] (a). No significant 
deviation of q from 1 is seen at 2-cr level. The best fit has in fact moved to Zq = 0.045 and 
q = 0.97, but is consistent with q been 1. This indicates that the putative Hubble bubble 
is perhaps an artifact of the data analysis. 

The CMAGIC method can shed more light on the monopole field. Not only because of 
its improved statistical errors, but also because it has a different sensitivity to systematic 
errors from extinction and light curve correction. The distance moduli of the sub-sample 
of SNIa derived in § [2] are fit with Equation (J3j). Here we did not apply a redshift cut to 
restrict the sample to only SNe at z > 2,500 km/sec as done in Ijha et al.l (120071 ) . Our 
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Fig. 5. — Correlations between the supernov a magnitu d e resi duals on the Hubble diagram 
and the host galaxy extinction as given by iJha et al.l (120071 ) . The linear fit results in a 
significant reduction of the x 2 /DoF and implies a correlation between the distance moduli 
and the extinction. Note that the linear correlation is determined mostly by the majority 
of the data with extinction below 0.5 mag, not those deviant points at larger than 1.0 mag. 
For clarity, the weighted average of the data points for Ay in the intervals [0, 0.1), [0.1, 0.2), 
[0.2, 0.3), [0.3, 0.4), [0.4 5), [0 .5,1.0), and [1.0, 2.0) are shown as solid circles. The open 
circles are from IJha et al.l (120071 ) . 



analyses show very little difference with or without this redshift cut. The results are shown 
in Figure El The conventional MLM yields contours that are very similar to the Ay corrected 
case shown in Figure H^b). In both MLM and CMAGIC, no monopole field more significant 
than 2-0" is detected. 

We thus conclude that our analyses do not support the presence of a monopole field 
in the local Universe. This brings the s upernova data into co nsistency with peculiar veloc- 
ity measured from clusters of galaxies (jGiovanelli et al.lll999l ). A similar result was found 
independently by Conley et al. (2007, in prep.) 
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Fig. 6. — Confidence levels of the monopole field from the conventional MLM method as 
given by equation (CD) (a, left panel), and from the CMAGIC method (b, right panel). The 
horizontal axes are the cut-off redshift, and the vertical axes the ratio of Hubble constant 
below and above z minus 1, as defined in equation [3j The contours show the 1-cr (solid 
lines), 2-0" (dotted lines), and 3-<r (dashed lines) significance levels. 



3.2. The Dipole Flow 



It is interesting to test wether there is a significant dipole field after the observed red- 
shifts are corrected to the reference frame at rest with the Cosmic Microwave Background. 
Such a dipole field could signify a global streaming motion of the local Universe, and ha s 
been found by previous studies (IRiess et al.lll995l ; IJha et al.l 120071 ; lHaugboelle et al.ll2006l ). 
The CMAGIC data, with the improved precision offer another probe of the local velocity 
field. 

Assuming that SNIa are perfect standard candles, the residual on the Hubble diagram 
is converted to peculiar velocity V pec using the formula V = zc(l — 10 a2<SM ), where z is the 
redshift of the supernova, and 5M is the residual magnitude on the best fit Hubble law. In 
Figure [7] we show the peculiar velocity field derived with this equation on the Aitoff diagram. 

Visual inspections show overall clustering of positive and negative values. The lumpiness 
of the velocity field is obvious in both Figure [3(b) for the CMAGIC Hubble diagram and 
Figure [7(a) for the MLM. Since none of the methods actually employs the position of the 
SN on the sky as input parameters, any angular correlation of the residuals are likely to be 
associated with unknown physical processes or systematic errors that are not yet accounted 
for in the analysis. 
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Fig. 7. — The peculiar velocity in the CMB frame on the Aitoff projection for the con- 
ventional MLM (a, upper panel) and CMAGIC (b, lower panel) method. The Galactic 
coordinate system is used. The red and blue symbols show positive and negative peculiar 
velocities, respectively. The size of the symbols are linearly proportional to the velocity and 
the thickness of the circles is inversely proportional to the velocity errors. The thick solid 
lines show the location of the Galactic plane. The inserts illustrate typical velocities in the 
figure. 
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Fig. 8. — The peculiar velocity on the supergalactic X-Y plane for the conventional MLM 
(a, left panel) and CMAGIC (b, right panel) methods. The star marks the location of the 
Great Attractor. The circles show positive peculiar velocities and the squares negative. The 
size of the symbols are proportional to the amplitude of the peculiar velocities. Note that 
in (a), the errors are large for the SNe at z > 0.025. 

The projected distribution on the supergalactic X-Y plane, shown in Figure [8] provides 
another view of the clustering of positive and negative peculiar velocity points. Figure [8] 
marks also the location of the hypothesized Great Attractor (GA). The supernova data do 
not reveal a convergence of the velocity field toward the GA. The results from CMAGIC, 
shown in Figure [8](b) seems to indicate excessive expansion velocity reaching far beyond 
the redshift of the GA. The positive and negative velocities appears to be distributed in 
two orthogonal directions. This pattern is less clear in Figure IE(a) for the conventional 
MLM method, due apparently to its the higher level of noise than the CMAGIC method. 
Nonetheless, as we shall see that the two methods actually reveal consistent information 
about the local peculiar velocity field. 

The data sparsely cover the entire sky, but with huge gaps that complicates spherical 
harmonic decomposition of the velocity field. To constrain the model parameters to a mean- 
ingful level, we impose further a constraint that the velocity field is axially symmetric. This 
is a rather arbitrary constraint but it reduces considerably the mathematical uncertainties 
and enables quantitative descriptions of the velocity field with simple mathematical formu- 
lae. We will consider both the dipole and quadrupole terms in this study. This is a natural 
step above the monopole field studied in the previous section. 
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For the axially symmetric dipole model, we model the velocity field in terms of the 
spherical harmonics Y®, 



V pec (l,b) = po + pifi = po + pi(cos&o cos 6cos(7o — I) + sin 6 sin 6), 



(5) 



where (7, b) is the longitude and latitude of the supernovae, p and p\ are the fitting coef- 
ficients of the monopole and quadruple field, (Iq, b Q ) are the coordinates of the symmetry 
axis, and fi = sin(b') with b' being the latitude measured in the polar coordinates with the 
polar axis pointing at (Zq, bo), i. e. the symmetry axis in Galactic coordinates. 



The x 2 °f th e model is then 



X 



Y,{Vi-V pec {liA)f/{8Vtl 



(6) 



where the errors SV of the peculiar velocity measurements are derived from the distance 
measurement errors. There are three major contributions to the velocity uncertainties: the 
measurement error, the error due to an unknown random velocity of the peculiar velocity 
field, and the error due to the intrinsic magnitude dispersion of supernovae. 

Applying equations (jSJ) and (jUJ) to the observed data, we found that a dipole fit re- 
sults in a reduction of \ 2 for both the MLM and the CMAGIC method. The resulting 
dipole for the entire sample is found to be at an amplitude of 293±100 km/sec and 289±94 
km/sec for the MLM and CMAGIC respectively. The direction of the dipole field is at 
(l ,b ) = (288°.0±^, -5°.7^ 6 7 2 ) and (Z ,& ) = (306°.0i^j, 17°.5^;|), for the MLM and 
CMAGIC, respectiv ely. These results are consistent with each other, and are also consistent 
with earlier studies (IRiess et al.lll995t IJha et al.l 120071 ; lHaugboelle et al.l 120061 ). This success 
is encouraging, and motivates further explorations of the peculiar velocity field from the 
CMAGIC method. 

A question of particular importance is the depth of the dipole pattern in the redshift 
space. If the Universe is homogeneous and isotropic on average, we should see a decreasing 
dipole amplitude as we average over larger and larger volume. To quantitatively investigate 
this problem, we have divided the supernova sample into three different sub-samples A, B, 
and C according to their redshifts. Each sample has roughly the same number of SNe. The 
samples A, B, and C are for SNe at redshifts below 0.01, between 0.01 and 0.025, and above 
0.025, respectively. The weighted mean redshifts (standard deviations) are 0.00675 (0.00239), 
0.01653 (0.00378), and 0.03450 (0.00887), for samples A, B, and C respectively. The optimal 
estimates of the parameters of these fits are given in Table [31 and are shown in Figure [9] and 
Figure [10J The dipole directions for the lowest redshift sample A is remarkably consistent 
with the direction of Norma cluster at (1, b) = (325°.3308, — 7°.2557 ) - the suspected center of 
the Great Attractor (see for example, iRadburn-Smith et al.l (120061 ) for a recent study of the 
GA). At higher redshift with sample B, neither the MLM data nor the CMAGIC data show a 
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prominent dipole. At even higher redshift sample C, the MLM data reveal no dipole, but the 
CMAGIC method does show a rather strong dipole whose direction is at (1, b) = (342°. 0^ 5 10 , 
11° S +21 1 



at an amplitude of 452^287 km/sec (Figure [TU1) . This is only a 2-cr detection, 



but it is a very intriguing possibility considering the fact that its direction is consistent with 
the dipole direction of sample A, and points right at the center of S hapley supercluste r of 
galaxies located at (1, b) = (327°, 32°) and a redshift of about 0.048 JSmith et al.ll2004h . 



Table 3. Dipole and Quadruple Fits to the Observed Peculiar Velocity Field* 
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*The errors are generally non-Gaussian, check Figures l9l 1101 [TT1 and 1121 for the error contours. 

^The superscripts M and c in the Model column stand for fits to MLM and CMAGIC data, respectively. The 
subscripts a, b> and c stand for sample A, B, C, respectively. The Model D, DQ, and Bip refer to the streaming 
motions as given in Equations \E\ [7] and \E\ respectively. 

^The redshift given in column < z > are the weighted average redshifts for the SNe in each sample. 
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Fig. 9. — The error contours of the parameters for the dipole fits using the MLM of SNIa. 
The left column shows the axis of symmetry (l ,bo), and the right column shows the fitted 
coefficients of the monopole and dipole components. From the top to the bottom, the panels 
show the fits to the entire sample, sample A, sample B, and sample C. For the contours of 
(l , bo), the probability distributions are marginalized over p and pi, for those of p and pi, 
they are marginalized over l and b Q . The x7 DoF 

were normalized to unity to make these 
contours. The solid, dotted, and dashed lines are 1-, 2-, and 3-cr confidence levels. 
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Fig. 10.— Same as Figure El but for the CMAGIC method. 
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It is also remarkable that the fits to the whole sample and samples A, B, and C are 
all consistent with the monopole component po being zero. This is a direct confirmation of 
the previous section where we found no measurable monopole components with a different 
approach. 

For sample A, the dipole fits reduce the \ 2 from 24 to 12 and 19 to 11, for the MLM 
and CMAGIC, respectively. Such a significant improvement of the goodness of fit suggest 
that the hypothesized model is indeed valid. After the dipole correction, the x 2 /DoF are 
well below unity for both MLM and CMAGIC. This can be understood if we have assumed 
a thermal peculiar velocity component that is too large. To make the x 2 /DoF to unity, the 
thermal peculiar velocity should be around 270 km/sec, assuming no intrinsic dispersions 
due to supernova standardization. At these low redshifts, a thermal peculiar velocity of 270 
km/sec introduces a magnitude dispersion of about 0.27 magnitude, which is well above ~ 
0.08 - the upper lim i t of th e intrinsic dispersions of CMAGIC calibrated supernova distances 



(jWang et al.ll2003bl . 120061 ). It thus appears that 270 km/sec is a reasonable estimate of the 
thermal velocity of the peculiar velocity field within the volume of z ~ 0.00675. This estimate 
is only weakly dependent on the intrinsic magnitud e dispersions of SNIa - subtracting the 



intrinsic dispersion of 0.08 mag (jWang et al.l 120061 ) results in an estimate of the level of 
thermal component of ~ 258 km/sec. 

The absence of a consistent dipole field for sample B, however, does not imply the 
absence of a coherent streaming motion at larger scales. It is possible that the data are still 
noisy, or the assumed mathematical model is wrong and is incapable of fully capturing the 
large scale streaming motions. Although it is true that future wide area supernova surveys 
will help in determining the high redshift bound of the dipole field, and its extension to 
higher redshift, we will investigate further with the current data set the nature of the dipole 
velocity component in the next sections by introducing other formulations of the coherent 
streaming flow. 



3.3. The Quadruple Flow 

The data allow for studies of higher oder components of the peculiar velocity field. 
To restrict the number of free parameters to minimal level, we will again assume axial 
symmetry so that the dipole and quadruple field point to the same direction. There is no 
physical motivation to model the velocity field this way other than for its simplicity. With a 
large enough data set, it should be possible to determine the complete quadruple components 
without the assumption of spherical symmetry. As we learned from the previous sections, the 
monopole field is small and non-detected. We will assume in this section that the monopole 
field is zero at all redshifts. 
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Fig. 11. — The error contours of the parameters for the dipole plus quadruple fits using the 
MLM data. The left column shows the axis of symmetry (Z , b ), and the right column shows 
the fitted coefficients of the dipole and quadruple components. From the top to the bottom, 
the panels show the fits for the entire sample, sample A, sample B, and sample C. For the 
contours of (/ ,&o)> the probability distributions are marginalized over pi and p 2 ; for those 
of pi and P2, they are marginalized over / and bo. The solid, dotted, and dashed lines are 
1-, 2-, and 3-a confidence levels. 
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Same as Figure El but for the CMAGIC method. 



-30 - 



In the reference frame where the polar direction is the symmetry axis, the model of the 
velocity field is then 

Vpec = pisinfr' +p 2 (3sin 2 6' - l)/2 (7) 

where b' is the latitude in the rotated reference frame and has a value in the range from —n/2 
to 7r/2. Just like the dipole model, the dipole plus quadruple model given in Equation (jjj) 
has also four free parameters: pi, p 2 , and the direction of the symmetry axis (l , b ). 

The results of the fits are shown in Figures [TTJ [T5J and in Table [3J As the figures show, 
the apex direction is at (l , b ) ~ (312°, 24°) for the full sample using the CMAGIC method. 
The MLM method gives results that are consistent with the CMAGIC method for the full 
sample, but with larger errors. We found also that for Sample A, the flow can be equally 
well described by a pure dipole (p 2 = 0) or a pure quadruple (pi = 0) with the CMAGIC 
method, but the significances of both the dipole and the quadruple field are less than 2 — a 
for both MLM and CMAGIC. For the intermediate redshift sample B, CMAGIC and MLM 
again give results that are in remarkable agreement; both methods show appreciable x 2 
improvement as compared to the pure dipole model. For sample C, only the CMAGIC 
method is able to reveal the smooth component of the peculiar velocity field. At less than 
2 — cr level, the flow of sample C shares the same apex direction and flow velocity as of 
sample A and B, although the best fit models are significantly different. The sample C of 
MLM method is unable to provide constraints on the flow field due to the larger associated 
errors. 



3.4. A Bipolar Flow Field 

It is found from §3.21 and §3.31 that there exists a coherent smooth streaming flow which 
requires a quadruple component. The exact functional form of this flow is unknown and the 
parametrization of the flow can be qualitatively different based on different model assump- 
tions. For example, at z below 0.01, a quadruple can provide an equally good fit as a dipole, 
but with the apex direction different. This is obviously a severe problem imposed by the 
lack of SNe in the zone of avoidance. Given this uncertainty, we explored further another 
functional form 

V pec = Pl cos[2(b' + p 2 )], (8) 

where the parameter p 2 is a phase offset such that when it is equal to zero the above equation 
gives a pure quadruple flow. The above equation reduces to a pure quadruple field when 

P2 = 0. 

The marginalized error contours of the fits are shown in Figure [13] and [HI for MLM 
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Fig. 13. — The error contours of the parameters for Equation [8] using the MLM of SNIa. 
The left column shows the axis of symmetry (l Q , fo ), and the right column shows the fitted 
coefficients of the velocity p% and the phase offset p 2 - From the top to the bottom, the panels 
show the fits to the entire sample, sample A, sample B, and sample C. For the contours of 
(/ , b ), the probability distributions are marginalized over p 1 and p 2 ; for those of pi and p 2 , 
they are marginalized over l and 60. The x7 DoF 

were normalized to unity to make these 
contours. The solid, dotted, and dashed lines show 1-, 2-, and 3-cr confidence levels. 
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and CMAGIC, repsectively. The optimal fit to the CMAGIC distances gives 

V pec = -340.0{^:f}cos[2(6 / + 24°.l{i|:|)}] km/sec (9) 

for the full sample, where the numbers inside the curly parentheses are the errors of the 
preceding numbers. To within 2-cr, the above equation, deduced from the full sample, applies 
also to all the sub-samples A, B, and C and to both MLM and CMAGIC. By removing this 
velocity field, the reduction of x 2 is from 90 to 62 for the CMAGIC Hubble diagram. The 
same operation with MLM results in a reduction in \ 2 from 78 to 63. Both fits are statistically 
significant. 

These model fits can be better illustrated in figures with the peculiar velocity plotted 
versus the latitude in the reference frame with the polar direction rotated to (Zo,&o)- To a 
good approximation, we can choose the direction (Z , &o) = (312°. 0, 25°. 7) to be the same 
for all the sub-samples. The fitting can be particularly simple as the a linear fit can be 
employed if the dependent variable is taken to be cos[2(6' + p 2 )]. The results are shown in 
Figures [TS] and [HJ Here we have again allowed for a monopole component, so the fitting 
formula is V pec = cq + c\ cos[2(6' + P2)] • The fitted parameters are shown in Figures [151 and 
[TBI and are tabulated in Table HI The MLM method and CMAGIC give overall consistent 
results, but the CMAGIC method shows much stronger signal especially at z above 0.025. 
These figures illustrate the need for more SNe at redshift above 0.025 to determine precisely 
the smooth peculiar velocity component and its redshift space extension. In the rotated 
reference frame, the peculiar velocity can also be plotted as a function of cos[2(6' + p 2 )] 
(see Figure [16]). A linear relation is expected if Equation (jSJ) is an acceptable description 
of the velocity field. This is indeed the case as shown in Figure [16] Figure [15] (c) and (d) 
shows the peculiar velocities from the CMAGIC method for SNe above and below z = 0.025, 
respectively. These figures illustrate again that the local streaming motion extends very deep 
in redshift space, to at least z = 0.025 (105 Mpc/hn) and, to a lower significance level, even 
the whole data set. 

Sample C has only one supernovae at b' below —30°. It would be very important to fill 
up this region. However, it is remarkable that the data points do agree well with an expansion 
law that describes reasonably well sample A and B (see Figure [ToT d) and Figure [TBT d)). This 
is an extraordinary property as the distance scales covered by these SNe is nearly 500 Mpc 
in diameter. 
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Fig. 14. — Same as Figure [T3J but for the CMAGIC method and equation [HJ It is seen 
from these figures that the peculiar velocity fields for samples A, B, and C share the same 
symmetry axis. 
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Fig. 15. — The residuals plotted against b', the latitude measured with the polar direction at 
(Z 0j bo) = (312°.0, 25°.7). The smooth solid lines show the fit V pec = 312.0 cos[2(&' + 24°.l)] 
km/sec, as deduced from the \ 2 fit to the peculiar velocity field of the full sample using 
the CMAGIC method, (a, upper left panel) shows the residual for the full sample of the 
conventional MLM method, (b, upper right panel) shows the same figure for the CMAGIC 
method. The residuals for the CMAGIC method using SNe below z = 0.025 and above 
z = 0.025 are shown in (c, lower left) and (d, lower right), respectively. No significant 
differences are identified in these fits, suggesting that the quadruple field extends at least 
beyond z = 0.025, possibly to z = 0.05. The hatted errors on the data points include only the 
observational errors, whereas the unhatted error bars include a 350 km/sec thermal peculiar 
velocity. 
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Fig. 16. — The same as Figure [T5| but the horizontal axis is now cos[2(6' +p 2 )]- The data 

points are expected to follow a straight line with the adopted model. The upwards triangles 

are for data with b' > — 24°. 1, the downward triangles are for data with b' < — 24°. 1. 
Table 4. Fixed Axis Axially Symmetric Fits to the Peculiar Velocity Field 



z range 




MLM 








CMAGIC 






DoF 


co (ct)* 


ci (<t)* 


Xo^ 


x 2 * 


co (<t)* 


ci (°0* 






- 0.08 


13.1 (50.4) 


-286.5 (76.2 ) 


78.7 


64.6 


41.4 (46.3) 


-377.4 (70.8) 


90.7 


62.1 


71 


- 0.025 


40.8 (55.1) 


-326.1 (84.2) 


49.3 


34.3 


62.2 (53.4) 


-347.2 (83.7) 


52.0 


34.5 


15 


0.025 - 0.08 


-05.4 (177.5) 


-145.5 (253.4) 


18.1 


17.6 


19.4 (91.0) 


-407 (130.9) 


36.2 


26.5 


21 



Xo an d X are X 2 before and after the fits. 



*In units of km/sec. 
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From Table HI we again did not find any significant monopole component at redshift 
above and below 0.05, when the symmetry axis is fixed to (Zq, &o) = (312°. 0, 25°. 7). 



3.5. Monte-Carlo Test of the Significance of the Streaming Velocity Field 

The uneven distribution of SNe on the sky is a worry of systematic errors that may lead 
to a false streaming velocity pattern. Extensive Monte-Carlo test is employed to assess the 
robustness of the streaming velocity field. We use the reduction of the x 2 as a measure of the 
goodness of the streaming pattern. In each Monte-Carlo realization, we keep the redshift, 
and (I, b) coordinates of each SN the same as were observed, but randomly assign a new V pec 
drawn from the observed peculiar velocities, with no repetition. The Monte-Carlo sample 
is fitted by the assumed streaming velocity models (Equations ([5]), (JTj), and ([8])) and the 
reduction of x 2 is recorded. Figure [17] shows the results for the streaming motion model 
given in equation 

From Figure [UJ we can see that the streaming pattern is significant at larger than 99% 
level for the full sample and for sample C, but is moderately significant at above 1 to 2-a 
level for the individual samples A and B. 

Note that the Monte-Carlo test were done independently for different sub-samples, so 
the joint probability of having a streaming motion pointing in the same direction by the 
various Monte-Carlo sample is significantly lower. This suggests that the observed streaming 
velocity field is highly significant. 

To test the effect of the uneven sampling on the errors of the parameter estimates, 
we have also applied bootstrapping simulations to study the parameter distribution. Only 
parameter space around the optimal fits given in Table [3] were searched for in the Monte- 
Carlo simulation. In this approach, the data are resampled 2000 times by redrawing SNe 
randomly from the original data set, allowing for repetition. Each sample is analyzed in 
the same way as the original observations. This generates 2000 Monte-Carlo realizations 
of the original observations. We have carried through the error terms in the bootstrap 
simulation; thus for each Monte-Carlo realization the resulting model fits have different 
errors associated with them. To calculate the bootstrap confidence level, we have weighted 
the fitting parameters by the inverse-squared of their corresponding errors to construct their 
probability distribution. The resulting constraints on model parameters are found to be in 
broad agreement with those from the x 2 fits shown in Table [3j The bootstrap test suggests 
that the smooth streaming pattern is significant at more than 3-cx level. The monopole term 
is consiste nt with zero f r om th e se analyse s . Thi s is in contradiction to the monopole field 



propose m 



Zehavi et ajj (119981 ); IJha et al.l (120071 ). but is consistent with the conclusions of 
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the previous sections of this paper. 



3.6. Systematic Errors That May Lead to the Observed Pattern 

The results derived here, although statistically significant, may subject to many sys- 
tematic errors. First, the data were acquired by different groups on different telescopes 
which studies supernovae only in certain part of the sky. The observations from different 
groups may have systematic errors as large as a few percent. Even for data taken by the 
same group with the same instrument, seasonality may introduce patterns across the sky 
as different parts of the sky are are accessible to the observer only in different seasons. To 
assess these effects, we studied the magnitude residuals as a function of the different research 
groups that acquired the original data. No dependence on the observers were identified. We 
conclude that the systematic errors introduced by different groups is unlikely the course of 
the streaming field, although it is difficult to rule out this possibility. 

Another source of systematic error is seasonality of the observations. Although never 
being reported to be the case, data taken in the winter may be systematically different from 
those taken in the summer, due to atmospheric changes or due to temperature changes that 
may affect the performance of the observations. However, no obvious pattern is observed, 
suggesting that this again is an unlikely origin for the large scale streaming field. 

We should stress that the streaming field is not to be due to the analysis method we 
used to homogenizing the supernova distances. Neither the MLM or CMAGIC method used 
in this study has a sensitivity to the position of the supernovae on the sky, and thus can not 
by itself introduce systematic correlations of the magnitude residuals on the sky. 

Sampling errors are certainly present and important. The presence of the zone of avoid- 
ance makes it difficult to observe supernovae close to the galactic plane. This causes am- 
biguities in decomposing the clustering patterns of the magnitude residuals on the sky into 
spherical harmonics. The dipole field for sample A, for example, is found to be pointing to- 
ward the suspected great attractor located behind the zone of avoidance. It can, however, be 
equally well modeled by a quadruple field pointing toward the center of Shapley supercluster 
(see Figure fT2l). More SNe are needed to disentangle this ambiguity. 

Systematic errors may also arise from corrections for Galactic dust extinction. The 
observed pattern may be the result of a similar systematic pattern in the dust extinction 
maps. The amplitude of the systematic error of the extinction map needs to be around 0.02 
mag to explain the amplitude of the streaming field for the sub-sample z at z > 0.025. 
However, if systematic errors of Galactic extinction is the cause, we expect the amplitude 
of the deduced coherent flow to be inversely proportional to the redshift of the supernova. 
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The analyses here for SNe below and above z = 0.025 give quite consistent expansion 
velocities at different redshifts. It thus seems that errors in Galactic extinction is not a 
favored explanation of the large scale streaming field. For the same reason, the observed 
velocity field is not due to an angular dependent Hubble constant. 



4. The Streaming Motion Corrected Hubble Diagram 

We can now apply the correction of the streaming motion pattern and construct a 
new Hubble diagram. We will use Equation ([9]) for all of the SNe, with the coordinate 
transformation given as 

b' = sin _1 [cos6 cos6cos(7 — I) + sin6 sin 6], (10) 

where (Iq, bo) = (312°. 0, 25.7°). After correction of the coherent flow given in equation (|UJ), 
the new Hubble diagram is shown in Figure [T9j As expected, the streaming motion correction 
dramatically improves the fits to the Hubble flow. Both the x 2 /DoF and the magnitude 
dispersions around the best fit Hubble line are significantly reduced. Most noticeably, the 
X 2 of the CMAGIC Hubble diagram decreased from 90 to 63, for 69 degrees of freedom, this 
makes the CMAGIC distances nicely agree with the a smooth Hubble flow. The magnitude 
dispersions on the Hubble diagram becomes 0.12 mag, the smaller ever for a sample of SNIa 
of this size. The correction of the streaming motion improves also the MLM Hubble diagram, 
at a more moderate level. 
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Fig. 17. — Monte-Carlo simulations of the significance of the streaming velocity field. The 
figures show the histogram normalized to peak value of 1 of x 2 reductions from the simulation. 
The solid lines show the results for CMAGIC, and the dotted lines for MLM. The vertical 
straight lines show the \ 2 reduction by fitting a streaming pattern as given in equation [8] 
to the MLM data (dotted) and CMAGIC data (solid). From top left clockwise, the panels 
show the results for the entire sample, sample A, B, and C, respectively. A total of 1000 
Monte-Carlo runs were performed for the full sample and each sub-sample. 
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Fig. 18. — Monte-Carlo bootstrap simulation showing the error levels of the parameters for 
equation [7J The numbers on top of each panel gives the mean and standard deviation of the 
corresponding parameter, for the full sample of the CMAGIC data. The errors are in good 
agreement with the results tabulated in Table EJ and in Figure [TU From top-left clockwise, 
the figures show the errors of pi, l , b , and p 2 , respectively. 
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Fig. 19. — Hubble diagrams after corrections of the coherent flow for the conventional method 
(a, upper panel) and the CMAGIC (b, lower panel). The Hubble diagrams shown here also 
have smaller magnitude dispersions as compared to those in Figure [31 A dramatic decrease 
of x 2 is found for the CMAGIC Hubble diagram as compared to Figure [31 
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Fig. 20. — The residuals of the coherent flow corrected distance moduli after the subtraction 
of the smooth Hubble flow as a function of redshift. The light dots show the residuals of 
individual SNe. The heavy dots show the average of data points in redshift bins (from left 
to right): to 0.01, 0.01 to 0.02, 0.02 to 0.03, 0.03 to 0.05, and 0.05 to 0.08. The solid lines 
show the Hubble line of a flat universe with Qm — 0.27, the dashed lines with Qm — 0, and 
the dash-dotted lines with Qm — 1-0. The panels are for the MLM (a, left) and CMAGIC 
(b, right) methods. 

To see the redshift dependence of the streaming field, we have plotted the magnitude 
residuals versus observed redshifts (see Figure [20] for the streaming velocity corrected Hubble 
diagram. The individual SN was grouped into smaller redshift bins [0, 0.01), [0.01, 0.02), 
[0.02, 0.03), [0.03, 0.05), and [0.05, 0.08). The weighted averages of the SNe in the different 
redshift bins are shown also in Figure [201 With the assumed thermal velocity of 350 km/sec, 
Figure [20] shows the data at different redshift are consistent with the assumed Q m = 0.3 
and f^A = 0.7 cosmology. An overall parabolic distribution of the average residual is 
indeed apparent in the data. We will defer the significance of such a pattern to future 
studies involving larger data set, but want to remark that at least in the redshift range 
below 0.01, the 350 km/sec is an over estimate of the thermal peculiar velocity as we found 
in $3] it is highly likely that a redshift dependent streaming motion will emerge from nearby 
SNIa observations. Two different cosmology models are also over-plotted Figure [233 which 
illustrates that it is potentially possible to study the nature of the dark matter and dark 
energy content of the Universe with SNe at redshifts as low as below 0.1. 
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Discussions 



The primary goal of this study is to reveal the large scale streaming motion of the local 
Universe. Earlier efforts based on distances derived from nearby galaxies and clusters of 
galax ies have shown results i n rough agreement with the current analysis. At redshift below 
0.02, iGiovanelli etaJ Jl998h found a dipole flow (in the CMB frame) of 200±65 km/sec 
toward (1, b) = (295°, 25°) from an I-band Tully-Fisher survey. iDekel et al.l (119991 ) found a 
velocity 370±110k m/sec toward (1, b) = (30 6°, 10°) from the Mark III velocity compilation. 
At higher redshift, lLauer fc Postman! (119941 ) found a 689±178 km/sec bulk motion (toward 
(1, b) = (343°, 52°) for an all sky survey of 119 clusters of galaxies to z ~ 0.03. iHudson et al. 



(119991 ) reported results from the Streaming Motions of Abell Clusters (SMAC) project, and 
found a bulk flow of amplitu de 630±200 km/sec toward (1, b) = (260±15°, -1±12°). More 
recently, IHudson et al.l (120041 ) analyzed 56 clusters within 120 h^ 1 Mpc in the SMAC sample 
and identified a bulk flow of 687 ± 203 km/sec toward (1, b) = (260° ± 13°, 0° ±11°), and 
found that flow does not drop off significantly with depth. 

The amplitude of these velocity fie l ds are cons istent with what w e found in this study. 
The apex direction of IGiovanelli et al.l (119981 ) and IDekel et al.l (119991 ) agree with the apex 
deduced from this study with Equation (JHJ). Due to the poor sky coverage of supernovae in 
sample C, the apex from the supernova data are quite sensitive to the underlying models. 
These differences may not be alarming considering the fact that so far most stu dies of large 
scale streaming motions have assumed dipole flows (cf. lHaugboelle et al.l (120061 ) ) . There are 
no strong physical reasons for using only a dipole field for the coherent flow, especially at 
very large scales. It would be interesting to see if fits using different model functions would 
bring the results derived from galaxy and cluster observations into better agreement with 
the supernova observations. 

The fact that the velocity field we derived points to the location of Shapley concentration 
at all of the redshift bins is a direct confirmation that Shapley supercluster contributes sig- 
nificantly to the streaming motion of the local Universe. The detection of a high-amplitude 
coherent flow on such a large scale argues for excess mass density fluctuation power at wave- 
lengths about 60-100/h Mp c, relative to the predictions of currently popular cosmological 
models (IHudson et al.ll2004r ). 



There may be a physical relation bet ween the quadruple field revealed here an d the 
quadruple fields identified from CMB maps (ITegmark et al.ll2003l ; iBennett et al.ll2003l ). The 
directions of CMB dipole, quadruple, and octopole are found to be aligned roughly in the 
same direction. Very likely these phenomena can be understood in a coherent picture. For 
example, it has been suggested tha t the CMB dipole may be caused by interaction of CMB 
photons with the local structures (IRakic et al.ll2006l ; llnoue fc Silk! 120061 ) . A comprehensive 
full sky nearby supernova survey is needed to fully address these questions observationally. 
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Such a survey can be performed on telescopes of large field of view, but not necessarily 
large telescopes. Future projects such as PANSTARRS and LSST provides some interesting 
scientific opportunity in this regard. However, for SNe below redshift of around 0.05, the 
supernova may be satu r ated w ith LSST. An array of small telescopes such as the system 
described in lYork et al.l (120061 ) is a powerful tool for mapping the local universe at z below 
0.05. 

An important question to ask is the level of the intrinsic magnitude dispersions of 
SNIa afte r standardiz a tion p rocedures are applied. This has been addressed in a previous 
paper by IWang et al.l (120061 ). We revisit this issue in light of the new understandings of 
the peculiar velocity field. We will consider only SNe at z > 0.01 to reduce the effect of 
the thermal component of the peculiar velocities. We then deduce the required amount of 
intrinsic dispersion in order to achieve a fit with reduced x 2 °f 1 f° r an assumed value of 
thermal component of the peculiar velocity. The conventional MLM method does not have 
sufficient statistical power to probe the intrinsic dispersion; for most reasonable values of the 
thermal component of the peculiar velocity, the amount of intrinsic dispersion needed is zero 
for MLM. For the CMAGIC method, this exercise yields values for the intrinsic dispersion of 
0.074 mag, 0.069 mag, 0.062 mag, 0.050 mag, and 0.029 mag for thermal peculiar velocities 
of 100, 150, 200, 250, and 300 km/sec, respectively. These are remarkably low values, 
especially when considering the fact that the data set assembled in this study came from 
diverse sources with unknown systematic errors. For the thermal velocity of about 250-300 
km/sec as we deduced from sample A, the intrinsic dispersion is found to be at the level 
of 0.03-0.05 magnitudes. This argues favorably that SNIa are indeed excellent standard 
candles, with the potential to reach distance precisions of around 2%. These numbers also 
suggest that the thermal velocities must be very low if the the intrinsic dispersions of the 
distance modulus deduced from CMAGIC is larger than 0.07 mag, i. e., the local streaming 
motion flow superimposed on the overall Hubble expansion, must be very smooth at the 
redshift range from to about 0.08. 

The reason that CMAGIC out performs MLM is perhaps related to the fact that 
CMAGIC weighs heavily on data about 2 weeks past optical maximum, whereas MLM uses 
data around optical maximum. At optical maximum, SNIa show more diverse behavior than 
at around day 15 past maximum. Spectroscopically peculiar SNIa such as SN 1991T are 
hardly distinguishable from normal SNIa at day 15 but the differences are obvious at around 
maximum. Furthermore, the polarization a r ound optical maximum decreases to practically 



unde tectable at day 15 (IWang et all Il996l . l2003al 120041 : iLeonard et al.l 120051 : IWang et al. 
20071 ). All these effects may have contributed. 



Finally, we want to remark the effect of the streaming motion on the determination of 
cosmo logical parameters using SNIa at z ~ 1. The large scale flow pattern implies that there 
could be correlated systematic errors at more than 0.02 magnitudes out to redshift of 0.1. 
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A sample of supernova evenly distributed on the sky is needed to average out the effect of 
the large scale correlated flow if no streaming motion correction is performed. 
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